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We perform a joint analysis of dwarf galaxy data from the Fermi Gamma-ray Space Telescope in 
search of dark matter annihilation into a gamma-ray line. We employ a novel statistical method 
that takes into account the spatial and spectral information of individual photon events from a 
sample of seven dwarf galaxies. Dwarf galaxies show no evidence of a gamma-ray line between 10 
GeV and 1 TeV. The subsequent upper limit on the annihilation cross section to a two-photon final 
state is 7 x 10 _26 cm 3 s _1 at 130 GeV, where the errors reflect the systematic uncertainty in 

the distribution of dark matter within the dwarf galaxies. 
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The search for dark matter annihilation directly into 
a photon final state is extremely important because the 
line emission occurs at an energy that corresponds to the 
mass of the dark matter particle (or thereabouts if the 
second particle is a heavy neutral particle) [IH7]. In ad- 
dition, line emission is free of background contamination 
as no known astrophysical process can result in line emis- 
sion at the energies of interest (a few GeV up to 10's of 
TeV). 

Recently there have been claims of the presence of a 
gamma-ray line at E 7 = 130 GeV [5rHU]. These studies, 
based on 3.5 years of data from Fermi, find a line emis- 
sion signature from the direction of the Galactic center. 
The interpretation of these results as dark matter anni- 
hilating directly to a photon final state implies a cross 
section of (av) s» [10~ 27 — 10~ 26 ] cm 3 /s. It is impor- 
tant to emphasize that this annihilation cross section is 
much larger than what one would expect from second or- 
der diagrams that lead to a two-photon final state (or a 
single photon and a Z gauge boson or h — for a sum- 
mary see e.g., [6]). Several dark matter interpretations 
for the alleged line feature have been offered [9J ITTlfrT] 
while other work raises doubts about the statistical sig- 
nificance of the line and its interpretation as dark matter 
[T51 fT§] . A recent search by the Fermi collaboration did 
not detect the presence of line emission in the Galactic 
halo (including the Galactic center) [20] . 

The Galactic center is clearly a place of interest when 
it comes to dark matter annihilation because of its large 
expected dark matter density [5T]. As the annihilation 
rate is proportional to the square of the number density 
of dark matter particles, its high density, coupled with 
its proximity to Earth, makes the Galactic center an at- 
tractive target for the search for an annihilation signal 
(e.g. [22]). 

In this short note we perform a search for dark matter 
annihilation to a photon final state in Milky Way dwarf 
galaxies using data from the Fermi Gamma-ray Space 
Telescope (Fermi). By virtue of their pristine dark mat- 
ter environment (absence of high-energy baryonic pro- 



cesses) and high concentration of dark matter, dwarf 
galaxies have been used to place the strongest bounds 
to-date on the s-wave annihilation cross section of dark 
matter [23l [21] • Given the paucity of background con- 
tamination along the lines of sight to the dwarf galaxies, 
it is natural to consider what limits the dwarfs may place 
on the annihilation cross section of dark matter into pho- 
ton final states. 

The approach we take is similar to [23l [25] . We per- 
form a line search by testing, at each line energy E 1 , the 
null hypothesis that the observed data was generated by 
background processes. Each hypothesis test is based on a 
test statistic T, which can be an arbitrary function of the 
data; however, it is vital that the choice of test statistic 
be made without reference to the data actually measured 
in the direction of the dwarf galaxies. 

We choose a simple form for the test statistic that com- 
bines the photon information from each of the dwarfs. 
Each photon % within a Region of Interest (ROI) of size 
1° is assigned a weight 10 based on which dwarf v it came 
from, its energy E, and its angular separation 9 from 
location of the dwarf. We denote this set of properties 
as Qi : {v, E,6}. The test statistic T is the sum of the 
weights of the photons detected within the ROIs centered 
on each dwarf: 

t = J2t„ (l) 

where the single-dwarf test statistic T v is 

T v = Y^w(Qi). (2) 

i=l 

Here, N v is the number of photons detected within the 
ROI centered on dwarf v. The weights w(Qi) and the 
total number of photons N v from each dwarf are random 
variables 1 . 



1 In previous work on the analysis of the continuum gamma-ray 
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To calculate the statistics of T it is useful to divide 
the parameter space of energy and angular separation 
(for each dwarf) into infinitesimal bins, each labeled by 
Q = {u, E, 6}. The number of photons detected in each 
bin is a random variable Xq. This total number of pho- 
tons is the sum of two random variables: the number 
of photons from dark matter annihilation Sq and the 
number originating from background processes Bq (i.e., 
Xq = Sq + Bq). It can be shown that the weight func- 
tion that maximizes an expected signal to noise ratio for 
a line emission search is |25j 

w(Q) = ^, (3) 

where sq and &q are the expected number of signal and 
background counts in the parameter space bin Q. 

Note that while sq and 6q are infinitesimal quanti- 
ties (being proportional to the size of the infinitesimal Q 
bin) their ratio is finite. In addition, sq depends on sig- 
nal characteristics in such a way that the expected signal 
in any infinitesimal Q bin is directly proportional to the 
annihilation cross section. Therefore, changing the anni- 
hilation cross section will simply scale the test statistic by 
a constant factor and will not affect any statistical con- 
clusions, i.e. this weight function is optimally powerful 
for any cross section. 

The quantity T v (Eq. [2]) is the sum of two terms: the 
weights of photons from dark matter plus the weights of 
background photons. These terms are independent vari- 
ables so the probability distribution function (PDF) for 
T v is the convolution of these individual PDFs. As in 
[23} [25] we model the background processes using data 
from the region surrounding each dwarf galaxy. The fun- 
damental assumption made is that the processes which 
give rise to the background nearby the dwarf also gener- 
ate the background at the location of the dwarf. 

For each dwarf we find the PDF of T v due only to back- 
ground processes by sampling the photons in the region 
within 15° of the dwarf. Sources from the second Fermi 
LAT source catalog [35] are masked with 0.8° masks (the 
95% containment angle for photons with energies greater 
than 10 GeV |27jV The sampling is performed by ran- 
domly placing 1° ROIs over the 15° field of view (reject- 
ing those ROIs which overlap with a masked source, the 
ROI centered on the dwarf, or the boundary of the field 
of view). The photons in these ROIs are then weighted 
according to Eq. [3] and summed as in Eqs.[l]and[2] 

In order to derive the PDF of T due to an annihilation 
signal, consider first a single dwarf T v as given in Eq. [2] 



emission from a combination of dwarf galaxies (23] . the weight of 
each photon was determined only by which dwarf it came from. 
The test statistic was therefore the weighted sum of the total 
number of photons collected from each dwarf. 



The quantity T v is the sum of N v independent, posi- 
tive random variables (the weights), where N v is drawn 
from a Poisson distribution with mean /Lt„, the expected 
number of dark matter photons from dwarf v. This dis- 
tribution is known as a compound Poisson distribution 
[25] , The PDF for T due to dark matter (Eq. [TJ for all 
the dwarfs is therefore the convolution of the individual 
compound Poisson distributions for each of the dwarfs. 
The PDF for each weight in the sum T v is the same and 
is found by dividing the energy-angular separation plane 
into infinitesimal bins and computing the probability that 
a detected dark matter photon will land in each bin. The 
weight assigned to a photon landing in each bin is set by 
Eq-i 

There are several methods for finding the PDF of the 
compound Poisson distribution T„ for dwarf v. An early 
algorithm was developed by Panjer [29] but we take ad- 
vantage of a straightforward and efficient fast Fourier 
transform (FFT) method [3D] which has also found use 
in astrophysics [3"TH3"3"] . 

For a single dwarf v, let T w „ be the Fourier transform 
(or characteristic function) of the probability distribution 
for the weight of a detected dark matter photon from v. 
The Fourier transform of the PDF for T v (due to dark 
matter annihilation), denoted Tt , is given by (see e.g. 

EDI). 

F Tv = cxp[^ (T WtL , - 1)]. (4) 

To incorporate both signal and background photons into 
the PDF for T we use the fact that a convolution is equiv- 
alent to multiplication in Fourier space. The full PDF for 
T is 

J T = [|cxp[/i J ,(J m ,„-l)] x (5) 

where J~b,v is the Fourier transform of the empirically 
measured distribution of the sum of weights due to back- 
ground processes for dwarf v. 

In practice the Fourier transforms are performed using 
an FFT on a discrete grid of possible T values. The 
single-event weight PDFs and the background PDFs are 
"tilted" [3U] before taking the FFTs to form J- W . v and 
Tb,v and the PDF of T is "tilted back" after applying 
the inverse FFT to Ft- The tilting prevents aliasing 
which can be induced by the FFT. 

The search for a line proceeds by first using Eqs. [5] 
& [3] to derive the PDF of T under the null hypothesis 
that there is no dark matter signal (all /x„'s are 0). The 
measured value of T, called T* , is obtained by summing 
the weights of all photons in the 1° ROIs centered on each 
of the dwarfs. The significance of the detection is the 
probability that T would be measured to be less than T* 
if the null hypothesis were true. For example, if there is 
99.7% chance that T <T* then a line has been detected 
at 99.7%, or 3a, significance. 



3 



The expected number of background counts bq is found 
by fitting a power law to all photons within 15° of the 
dwarf (excluding the central 1° and the masked sources). 
For purposes of weighting, the background is assumed 
to be statistically isotropic, i.e. independent of angular 
separation from the ROI center. This may not be true in 
practice due to the presence of unresolved sources; how- 
ever, the background sampling automatically includes 
any non-Poisson aspect of the background in the PDF 
of T v (or equivalently in Tb,u)- 

The expected number of dark matter annihilation 
events sq , detected from a particular dwarf, with energy 
between E and E + dE and with angular separation in a 
solid angle interval d£l(9) is 



J 



(<rv) dN 1 (E) 



8nM x 



dE 



e(E)PSF{E,9)dEdn(9). (6) 



In the above M x is the mass of the dark matter particle, 
(<jv) is the velocity-averaged annihilation cross section 
into a pair of gamma-rays, and dN^/dE is the number 
of photons per energy interval emitted per annihilation. 
The point spread function PSF(i?, 9) is the probability 
per solid angle of detecting a photon of energy E an an- 
gular distance 9 from the source, and e(E) is the detector 
exposure in units of cm 2 s. The quantity J quantifies the 
dark matter distribution within a particular dwarf (341 — 

SI]. 

We use the publicly available data from the Fermi 
Science Support Center (FSSC) 2 and version v9r27pl 
of the Fermi Science Tools. We extract all photons 
of evclass = 2 using the tool gtselect in the Mis- 
sion Elapsed Time interval [239557417 - 357485329] in 
the energy range between 8 GeV and 1 TeV, and with 
zmax = 100. We select good time intervals (with all 
standard recommendations as stated on the FSSC), 
and compute e(E) and PSF(I?, 9) using gtpsf with the 
P7S0URCE_V6 instrument response functions. 

The dark matter annihilation is modeled as point 
source emission from each dwarf, and we utilize the values 
for J given in |24j . For a line search the energy dispersion 
of the detector can be important. We incorporate this un- 
certainty by giving a width to dN^/dE. The spectrum 
due to line emission is simply dN 1 /dE = 25(E — E 7 ) (but 
see also 0H2)- The 68% containment on the energy un- 
certainty of Fermi is approximately 10% for photon ener- 
gies above 10 GeV. We model this energy uncertainty by 
setting the annihilation spectrum dN^/dE to be a Gaus- 
sian centered on M x , normalized to 2, with a standard 
deviation of 10% of the mean. We have reproduced the 
analysis with top hat distributions with widths from 5% 
to 30%. The effects are small and leave our conclusions 
unchanged. 
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FIG. 1: A snapshot of every photon having an energy between 
10 GeV and 1 TeV that has been detected within 1° of each of 
the seven dwarfs. The black dashed line is at 130 GeV [5] [9]. 
The gray region is ±15% around 130 GeV, a rough gauge of 
the energy dispersion of the LAT. The vertical axis is scaled 
according to solid angle so that an isotropic distribution of 
photons will be spread uniformly along this axis. 



Figure [T] shows the individual photon events between 
10 GeV and 1 TeV that were detected within 1° of each 
of the seven dwarfs. The vertical axis measures the an- 
gular separation between the event and the center of the 
dwarf. It is scaled according to solid angle so that an 
isotropic distribution of events should be distributed uni- 
formly over the vertical axis. There are no photons with 
energy within 15% of 130 GeV (gray shaded region). The 
68% energy resolution of the LAT ranges from about 8% 
at 10 GeV to about 14% at 1 TeV while the 68% con- 
tainment angle (PSF) varies from 0.3° to 0.2° over this 
energy range (95% containment is about 0.8°) [37]. We 
conclude that the dwarfs show no evidence of a gamma- 
ray line at 130 GeV. 

Using the formalism described above we perform a 
search for line emission over a range of energies. A Gaus- 
sian energy spectrum with a standard deviation of 10% 
is used to calculate sq (Eqs. [3] and [6|. We perform a 
separate search for each possible line energy, taking 100 
log-spaced steps from 10 GeV to 1 TeV. The results of the 
search are illustrated in Fig. [2] Note that the inclusion of 
a trials factor dilutes the significance of any line. We can 
make a very rough estimate of the number of "indepen- 
dent" trials by assuming that a search for a line at E 1 
uses the photons in the window E^(l±a). If the energy of 
the (—a) edge of the window is E\ the upper edge of the 
window is at an energy E 2 = Ex (1 + a)/(l — a). There- 
fore the number of "independent" (i.e. non-overlapping) 
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FIG. 2: Results of the a search for line emission using an opti- 
mized combined search of seven dwarf galaxies. The horizon- 
tal axis represents the energy of the gamma-ray line searched 
for. The left vertical axis is the significance of the detection 
(in terms of Gaussian standard deviations). The right vertical 
axis incorporates a trials factor of 24, roughly the number of 
independent energies searched. The non-significant peak at 
200 GeV is due to a single photon from Sculptor (see Fig.fl]). 
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FIG. 3: 95% upper limits on (av) for annihilation into a pair 
of photons each having energy £? 7 . The black line is the limit 
using the best fit J values for the dwarfs. The blue region cor- 
responds to the 95% systematic uncertainty in the estimates 
of J. The two points are the dark matter interpretations for 
the tentative signals observed by [51 [9] under the assumption 
of an Einasto dark matter profile and annihilation into two 
gamma-rays, with 95% error bars. 



specified by E max = E min [(l + a) /{I - a)] n . An energy 
window of a = 0.10 corresponds to about 24 trials. On 
the right vertical axis of Fig. [2] we plot the significance 
including a trials factor of 24 as a rough guide to the true 
significance of any tentative line. It is clear that the data 
do not strongly suggest that line emission is present at 
any energy. 

Given that there is no evidence of line emission from 
the dwarfs we can place upper limits on the annihilation 
cross section into two photons. In this case, the weight 
choice analogous to Eq. [3] that maximizes the signal to 
noise ratio is w(Q) — sq/(6q + sq) [25] . For each mass 
we find the cross section above which there is less than a 
5% chance of measuring the test statistic T to be smaller 
than observed. The resulting upper limits are plotted in 
Fig. [3] (together with the results from [H [H])- By far, 
the largest source of systematic uncertainty is in the J 
values for the dwarfs. The black line in the figure is the 
limit found when the J values are set to their best fit 
values found in [24] . The effect of varying the J values 
within their observational uncertainties is shown by the 
blue shaded region. One at a time, we set the J value 
for each dwarf to its upper or lower 95% error bar and 
recompute the 95% cross section upper limit. The dif- 
ferences induced by each dwarf are added in quadrature 
to produce the boundaries of the shaded region. This 
procedure gives an estimate of the systematic effect due 
to the difficulty of determining each dwarf's dark matter 



distribution. 

For annihilation channels producing continuum emis- 
sion (e.g. into heavy quark or lepton pairs) dwarf galax- 
ies provide strong limits on the annihilation cross sec- 
tion [23 Ell EH [391 SHU]. It is challenging to pro- 
duce such limits from the Galactic center: despite the 
high dark matter density (J value hundreds to thousands 
of times larger than the dwarfs) the astrophysical back- 
ground cannot be easily subtracted or modeled. However, 
a gamma-ray line search is not hindered by these back- 
grounds to the degree that a continuum search is. For 
this reason, the Galactic center may be a more attrac- 
tive target when searching for line emission. The upper 
limits obtained by E2] are much stronger than those 
obtained here from the dwarf data. A recent search by 
the Fermi collaboration for gamma-ray lines in the Galac- 
tic halo (including the Galactic center) [2U] did not show 
evidence for a 130 GeV line and places stronger upper 
limits than found here. 

It appears that the large increase in dark matter den- 
sity, and the proximity of the Galactic center are much 
more constraining than are dwarf galaxies when it comes 
to line emission searches. At the present time dwarf 
galaxies can neither confirm nor deny a dark matter line 
interpretation of the Galactic center data. 
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